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Abstract 

We construct a Schwarzschild-type exact external solution for a the- 
ory of gravity admitting local Galilean invariance. In order to realize 
the Galilean invariance we need to adopt a five-dimensional manifold. 
The solution for the gravitational field equations obeys a Birkhoff-like 
theorem. Three classic tests of general relativity are analyzed in detail: 
the perihelion shift of the planet Mercury, the deflection of light by the 
Sun, and the gravitational redshift of atomic spectral lines. The Galilean 
version of these tests exhibits an additional parameter b related to the 
fifth-coordinate. This constant b can be estimated by a comparison with 
observational data. We observe that the Galilean theory is able to repro- 
duce the results traditionally predicted by general relativity in the limit 
of negligible b. This shows that the tests are not specifically Lorentz 
invariant. 
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1 Introduction 



In recent papers, we have examined various aspects of Galilean covariance in a 
flat spacetime (for instance, see Ref. [T]), previously discussed in Refs. E]-!!]. 
This formalism allows us to express physical theories in a tensor form by em- 
ploying a 5-dimensional flat manifold characterized by a metric tensor (given in 
Eq. Therefore, as mentioned in these early papers, the reduction from this 

extended space can be done to the (3+1) Galilean spacetime or Minkowski space- 
time, thereby providing a unified framework to treat both kinematics. Among 
our motivations for the project of extending the study on Galilean symmetry to 
gravitation is the following statement (p. 92 of Ref. [S]) : "... general covari- 
ance does not imply Lorentz invariance - there are generally covariant theories 
of gravitation that allow the construction of inertial frames at any point in a 
gravitational field, but satisfy Galilean relativity rather than special relativity in 
these frames.'^ Our intent is to examine various consequences of replacing local 
Lorentz invariance with local Galilean invariance. Traditionally, an obstacle to 
the implementation of Galilean invariance in a four-dimensional manifold is the 
absence of a metric tensor satisfying this requirement locally [6] , but as we will 
discuss, a metric structure exists if we extend the spacetime to (4 -I- 1) dimen- 
sions. As far as we know, a systematic study of gravitational theories which 
admit Galilean invariance has never been performed before. 

Let us begin by briefly describing the Galilean-covariant formalism (more 
details can be found in Refs. The first step consists in constructing 

Lorentz-covariant action functionals defined on a fiat Galilean spacetime, defined 
as a (4 -t- l)-dimensional flat Minkowski manifold with light-cone coordinates, 
defined by Eq. ([3]). A Galilean five-vector, such as (x, x^, x^), transforms under 
Galilean boosts as follows: 



where v is the relative velocity. The first two lines of Eq. ([T]) define the well- 
known pure Galilean transformations in (3 + 1) dimensions. 

Note that a general Galilean transformation, in the (4+1) Galilean space- 
time, has a form similar to the Poincare transformation, x'^ = A'^+x" + a" . 
The tensor methods utilized in the Galilean-covariant formalism are based on 
the scalar product. A • _B = A • B — A^^B^ — Ar^B^, which is invariant under the 
transformations in Eq. ([T]) . This suggests the following Galilean metric in order 
to describe the flat Galilean spacetime: 
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This defines light-cone coordinates which, through the change of coordinates, 

_ _4 ~\~ X _^ X X . . 

can convert the metric, Eq. ([2]), into its diagonal form: 

7?^, =diag(l,l, 1,-1,1). (4) 

Therefore, dimensional reductions from (4 + 1) theories can lead to either a 
Lorentz-covariant or a Galilean-covariant theory in (3 + 1) dimensions [5]. 

The projection from the extended manifold onto the (3 + l)-dimensional 
spacetime can be defined as 

where x is the position vector, t is the time, c is a parameter with the dimensions 
of velocity, and S is the extra coordinate. The five-momentum 

Pp, = id^ = (iV, idt, ids) = (^P' ' (^) 

where p is the three-momentum, E the energy and m the mass, suggests that the 
additional coordinate, may be seen as being conjugate to the mass. Clearly, 
the theory in the (4 -I- l)-dimensional manifold is Lorentz-invariant since the 
two metrics are equivalent. It is through the definition of the physical quanti- 
ties that the reduction to a (3 -I- l)-dimensional spacetime leads to a Galilean 
theory. From a group-theoretical point of view, this dimensional reduction is 
motivated by the fact that the Poincare group in (4 -I- 1) dimensions contains, 
as a subgroup, the extended Galilei group of (3 -I- l)-dimensional spacetime. 
Indeed, the Poincare group of (4 + 1) spacetime is generated by 15 elements (10 
Lorentz transformations, M^^, and 5 translations, P^, where fi,!^ = 1, ... ,5). 
The 11 generators of the extended Gahlei group consist of the 3 rotations {M12, 
-^23, -^31), 3 Galilean boosts (M41, M42, M43), 4 spacetime translations {Pi,P2, 
P3, P4), and 1 central charge, P5, which is the mass. 

The Galilean symmetry is important because it is the kinematic symmetry 
which underlies non-relativistic, or low-energy, phenomena [7]. By simply tak- 
ing low-velocity approximations of relativistic, Poincare-invariant, expressions, 
we do not necessarily obtain theories with suitable symmetry properties. More- 
over, such approximations may not provide a complete description, as we will 
illustrate hereafter. Also, a careful consideration of Galilean invariance proves 
that some concepts or properties usually thought to be specifically relativistic 
can in fact be shared with non-relativistic systems (e.g. the spin of particles 
[S]). Hereafter, we will reach a similar conclusion with the classic tests of general 
relativity (GR). 

These points are illustrated clearly by the theory of Galilean electromag- 
netism [9l [1^. Firstly, the electric field E and magnetic field B transform. 
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under a Lorentz transformation, as follows: 

E' = 7(E + V X B) + (1 - 7)^, 

B'=7(B-lvxE) + (l-7)^, 

where the parameter v represents the relative velocity, and 7=1/ \/l — w^/c^ 
is the Lorentz factor. An inconsistent approximation (see remarks in Ref. [3]) 
is given by the relations: 

E' = E + vxB, (6) 
B' = B - X E. (7) 

For a recent and clear discussion, see Ref. [10]. Moreover, Le Bellac and Levy- 
Leblond found that there exists not one but two consistent Galilean limits of 
electromagnetism; the "electric" and the "magnetic" limits [3]. 

It is noteworthy to mention that Galilean invariance and the equations above 
arc not simply a theoretical figment, since they describe concrete physical sys- 
tems. Indeed, the two Galilean limits of the Maxwell's equations, given above, 
are the basis of electroquasistatics and magnetoquasistatics [TT]. In turn, they 
are connected to electrohydrodynamics and magnetohydrodynamics, which have 
many applications in geophysics, astrophysics, physics of fluids and plasmas, en- 
gineering, etc. 

The purpose of this paper is to construct an effective theory of gravity by 
performing a slight modification to the principle of general covariance that fol- 
lows from the equivalence principle. Quoting Weinberg once again (p. 91 of Ref. 
[S]): "/i [the Principle of General Covariance] states that a physical equation 
holds in a gravitational field, if two conditions are met: 1 . The equation holds 
in the absence of gravitation; that is it agrees with the laws of special relativity 
when the metric tensor g^p equals the Minkowski tensor rja^ and when the affine 
connection T^^ vanishes; 2. The equation is generally covariant; that is, it pre- 
serves its form under a general coordinate transformation x — *■ x'" . Henceforth, 
we propose to replace the condition 1 by, "1'. The equation holds in the absence 
of gravitation; that is, it agrees with the laws of Galilean relativity when the 
metric tensor gajs equals the Galilean tensor 770^3, Eq. ([2]), and when the affine 
connection F^^ vanishes." The condition 2 is kept without change. This new 
version of the principle of general covariance, with conditions 1' and 2, gives 
rise to a general (globally) covariant theory of gravity that is locally Galilean 
invariant. It is in this sense that the expression "Galilean covariant" should 
be understood henceforth. Wc shall refer to this new theory of gravity as the 
Galilean gravity. 

We call the Galilean gravity an "effective theory" in the sense that wc expect 
it to be valid in the same range of distances (in the tests described in Section [3]) 
as the post-Newtonian approximation (see Chap. 9 of Ref. [5]). Nevertheless, 
we will obtain an exact solution using local Galilean invariance (see Section 
[2]) rather than the expansion typically encountered with the post-Newtonian 
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approximation. It should be emphasized that, although we work in a (4 + 1)- 
dimensional manifold, our approach and purpose are entirely different from 
other 5-dimcnsional theories (e.g. Kaluza-Klein models, de Sitter universe, or 
the 'bulk' utilized in string theory). 

In 1916, Schwarzchild published a solution of Einstein equations for the 
static spherically symmetric gravitational field [12]- [2]. In Section [21 wc dis- 
cuss an analogous solution in the Galilean context, called 'Galilei-Schwarzschild 
solution'. This solution is exact, unlike other effective theories that ensue from 
approximations or expansions. Our solution contains the same event horizon 
as the familiar Schwarzschild metric. We also find the Galilean analogue of the 
Birkhoff theorem. In Section [31 predictions of this solution are confronted with 
three classic tests of GR. In Section 13. 1[ we find that the Galilei-Schwarzschild 
solution describes the perihelion shift of Mercury as does its counterpart of 
GR. In Section 13. 2( we observe that it predicts the defiection of light rays by 
massive bodies like our Sun. In Section 13. 3|, we obtain a similar result for the 
gravitational shift of atomic spectral lines. The results of these tests depend 
on an additional parameter 6, which is related to the fifth-coordinate. This 
constant b can be calculated by comparing with observational data. Wc find 
that the Galilean theory can reproduce the results predicted by GR in the limit 
6 ~ 0. Thus the tests are not specifically Lorentz invariant. Section [H contains 
concluding remarks. 

2 Galilei-Schwarzschild solution 

The principle of equivalence states that in the presence of a gravitational field, 
we can choose, at every spacetime point, a locally inertial coordinate system 
where the laws of nature take the same form as in a non-accelerated Cartesian 
coordinate system in the absence of gravitation. As mentioned in p. 91 of 
Ref. [S], this can be restated as the Principle of General Covariance, which 
requires the equations of Physics to preserve their form under general coordinate 
transformations, meaning that these laws are properly written in a tensorial 
form. 

As we build the 5D Galilean theory for gravity upon the same requirement of 
general covariance demanded by GR, it is only natural that the field equations 
for Qfj^i, in the context of our model are formally the same as in the usual GR. 
The manifold is equipped with a metric tensor g^^ which satisfies 




(8) 



where the Ricci tensor components, 




(9) 



(10) 
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2.1 Galilean metric in curved spacetime 

In GR, in order to derive the Schwarzschild solution, wc begin with a 4-dinicnsional 
metric tensor of the form [T51 



V 














sin^ I 













(11) 



It exhibits spherical symmetry in the spatial sector. In the relativistic theory, 
the functions A (r, t) and v (r, t) are calculated by (i) substituting this ansatz 
into Eqs. ([HI); (ii) considering T^i, = in Eq. ([S]), corresponding to the solution 
outside the massive source of the field; and (iii) imposing the boundary condition 
limr_»cx3 A = limr_,oo = 0, in order to recover the Minkowski space as a limiting 
case, since far from the source the spacetime should be flat. The Birkhoff 
theorem assures that A = A (r) and v = (?') (sec, e.g. Ref. [H]). 

Just as Eq. ^ converts the flat Galilean metric of Eq. ^ to the Minkowski 
metric given in Eq. ^ , the Galilean version of the Schwarzschild solution should 
be related by a similar change of coordinates to Eq. pT|) . with one trivial 
additional component. 
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(A more general ansatz exhibiting spherical symmetry would be 



(12) 
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(13) 



This reduces to Eq. ([T^ when we take fj, = 0. This choice is not only adequate 
because of its simplicity but also due to the fact that it leads to a Schwarzschild- 
type of metric. Other solutions could be obtained with a non-trivial /z.) The 
metric in Eq. leads to a Minkowski-like metric in (4+1) dimensions under 
the constraint linv^cxj A = limr_>oo 1^ = 0- The relativistic limit of a Galilean 
theory is obtained from a light-cone-type transformation in spherical coordi- 
nates, 

=4 j_ ;f.5 =4 =5 

r = f, d = d, ^ = (^, x^ = —^^ = (14) 

Therefore, in order to retrieve the general Galilean metric, we must transform 
Eq. p2[) to the unprimed reference frame through the usual transformation rule 



5 



for a tensor of rank two, 



and the metric in Eq. (fT^ becomes 



( 
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sin^ I 



dxP dx" _ 








(15) 



1) 
1) 



1) / 



Henceforth wc sUghtly gcnerahzc this metric as follows: 
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r sm I 
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r, X 



x^) 
x^) J 



(16) 



This must obey the following constraints, which ensure the rclativistic limit 
and follow from Eq. (I15p for (744 and 555: 

1 



h-2g) = l 



(17) 
(18) 



Let us justify this statement. In the limit r —^ 00, the gravitational field should 
become negligible, so that the metric tensor reduces to the flat Galilean metric, 
Eq. (12). The transformation, Eq. (|14p . allows us to obtain the rclativistic limit. 
Since Eq. (|14p is valid in the limit r 00, it can be applied everywhere. This 
assumption is consistent with the principle of minimal coupling, a prescription 
used in GR to express the laws of physics valid in flat spacetime in a covariant 
fashion. Under the transformation in Eq. p4|) . the metric becomes 

(5^^)4x4 
1 



where (gfn')^^^ is the standard Schwarzchild metric of GR. In this case, (i) the 
component 355 of the transformed metric, 355 = ^ (^44 — 2^45 + 355), should be 
equal to 1; and (ii) the component 345 = ^ (344 — ^55) would have to be null. 
These facts lead to Eqs. ((T7)) and (fT8)) for the metric, Eq. (fTB]) . Under these 
constraints it assumes the form: 



/ eA(.,.^.^) 



\ 



r 






sin^ I 



-l + /(r,x4,x5) 



-1 + f {r,x'^,x^) f{r,'. 
This is the ansatz employed in the following for the Galilean metric. 



(19) 
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2.2 Static spherical vacuum solution 

The field equations in vacuum, for which T^^ = in Eq. ([5]), reduce to 

R^. = 0. (20) 

In this section, wc analyze the vacuum field equations for which the arbitrary 
functions in the metric tensor do not depend on the a;^-coordinate; that is. 



X^X{r,x^) ; f = f{r,x^) 



(21) 



We call this metric a static solution, by analogy with the GR context where the 
coordinate is the time. Henceforth, we denote the derivatives as follows: 



f=drf; f=d5f. 



(22) 



We can calculate the explicit form of the Einstein field equations, Eq. (PO]) . 
for the metric in Eq. with the simplifying hypothesis, Eq. (|2ip . after ob- 
taining the inverse matrix g'^" , the corresponding Christoffel symbols, Eq. pO]) . 
and finally the Ricci tensor, Eq. ([9]). The resulting non-zero Ricci components 
are 



Rii = e 



^ 2 (1 - /) / A + (-/) (-1 + 2/) + 2a) 

4(1-2/)' 

V ^ 2^ + (-l + 2/)(rA'-2r) 
r 2(1-2/)' 

2/ /' + (-! + 2/) (Xf'-2f' 



4(1-2/)^ 



Ri 



_ 2//' + (-l + 2/)^A/'-2/'j ^ A 

-Ri5 — ; "2 \ — Rli 

4(1-2/)' r \r 



R22 — 1 



A' 



2/' 



r (1 - 2/) 



R33 R22 sin 



(23) 

(24) 

(25) 

(26) 
(27) 



R. 



44 



2/ /' - (-1 + 2/) / + (-/) (-1 + 2/) / A 



4(1-2/)^ 

2r' + (-i + 2/)(fv-2/") r 

4(1-2/) 



(28) 
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R45 — 



2 (-1 + /) P - (-1 + 2/) / + (1 - /) (-1 + 2/) / A 



4(1-2/)^ 
2/'2 + (-1 + 2/) (/'A' - 2/") , /' 



4(1-2/) 



+ 



(29) 



i?5 



2/ P - (-1 + 2/) / + (2 - /) (-1 + 2/) / A - (-1 + 2ff ( A2 + 2A 



4(1-2/) 

2/^ + (-i + 2/)(ry-2ro r' 

4(1-2/) 



(30) 



The condition in Eq. (j20[) requires that each component be equal to zero. 
Then. i?i4 = in Eq. ([M]) . By using this information in Eq. (I?51) . we find 



A = 0. 



(31) 



Therefore, the function A (r, x^) depends solely on r: A = A (r). 

From Eqs. (|24p and and the fact that A = 0, together with the require- 
ment that 1 — 2/ ^ 0, we find 



//'-(-l + 2/)/' = 0. 



(32) 



Also, since the Ricci components are equal to zero, we must have + R55 
2i?45 = 0, so that Eqs. ^ to ^ lead to 



/2-(-l + 2/)/ = 0. 



(33) 



It is important to observe that, with Eqs. ([?T|) . ([5^ and we find that 
Eqs. to ([50]) . for the x^- independent solution, reduce to 



P A' 2f^ + (-l + 2/)(fy-2/") 

illl — 1 :^ = U, 

2(1-2/)^ 



R22 — 1 



i?44 = ^45 = ^55 — — 6 



A'-2 



2/' 



= 0, 



r 1-2/ 
i?33 = i?22 sin2 e = 0, 
2/'2 + (-1 + 2/) (/'A' - 2/") , /' 



4(1- 2/) r _ 

With the condition 1 — 2/ 7^ 0, we can write Eq. ([57]) as follows: 

2f^ + (-l + 2/)(fy-2fO ^ 2f 

2(1-2/)' r(l-2/)- 

This can be substituted into Eq. ((34|) and it leads to the relation 

A' -0 



(34) 

(35) 
(36) 
0. (37) 



(38) 
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which is valid for r 7^ 0. 

If we take the derivative of Eq. ((38|) with respect to , and use Eq. ((3T|) , 
we obtain 

(1 - 2/) /' + 2/ /' = 0. (39) 

This relation, together with Eq. allows us to prove that / does not depend 

on . Indeed, multiplying Eq. p2|) by 2 and subtracting from ((39l) gives 

(-l + 2/)/' = 0. 

This implies that either 

/' = 0, (40) 

or 

-l + 2/ = 0. 

This last possibility is prevented by consistency of the theory since the factor 
— 1 + 2/ appears in the denominator of the Ricci components, Eqs. ([^5)) to ([50]) . 
From Eq. (pS)) . we find the relation 

A + In (1-2/) =lnD, 

where D depends on x^ only. We use In I? for convenience. Indeed, this equation 
relates A and / as follows: 



D 



e 



(41) 



1-2/ 

If we substitute Eqs. ([35]) and (gll) into Eq. we obtain 

(./)' = ^ ^ 

where £^ is a function of xJ' only. Therefore, 

rf = Er + F, 

with F depending on x^ only. Since = (1 — D) /2, we can write 

Up to this point, we have / — f{r,x^), D = D{x^) and F — F{x^^). By 
differentiating Eq. (|42p with respect to r and a;"", we find 

From Eqs. (|33]) and (gD]), we conclude that 

F 0, (44) 
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so that F is constant. Therefore, 

/ = -f- (45) 
By substitutmg Eq. (gO]) into Eq. ([32]), we find 

//' = 0. 

Since Eq. (|42|) prevents /' from being zero, this equation gives 

/ = 0. (46) 

Note that with Eq. the constraint in Eq. (|33p is autoniaticaUy satisfied. 

In addition, Eq. ((45|) gives 

I) = 0, (47) 

so that D is also a constant. 

According to Eqs. ((3T|) and ([46]), we see that 



A-A(r), / = /(r), 
thereby obtaining the following theorem: 



Galilei-Birkhoff theorem: A static spherically vacuum solution is necessarily -independent. 



This theorem is subjected to the restriction specified in Eq. p3)) : that is, 
our ansatz for the Galilean metric is the one consistent with Eq. (jl2[) . 

Since the metric of the fiat Galilean space-time, Eq. ([2|), is symmetric under 
the exchange of and , we obtain the usual Birkhoff theorem for the time 
x^; that is, if we begin with a solution independent of x^, then the solution is 
necessarily static. 

The metric is determined up to D and F, which are both constant, due to 
Eqs. ([37]) and (gj]), respectively. By inserting Eq. (gH) into Eq. (gT]), we obtain 

r 

The constant D is determined by the behavior of the metric at large distances 
from the source. In order to ensure that the metric reduces to the Galilean 
metric, Eq. in the limit r oo, we must impose lim,,^oo / ('') — 0. From 
Eq. (g5]), this leads to 

D = l. (49) 

The constant F is determined from the relativistic limit in the weak-field 
approximation. The light-cone transformation, Eq. (HH), and Eq. give 
944 = (544 + 2^45 -f- 355) /2. For the Galilean metric, Eq. (|16p . subject to the 
constraints in Eqs. p7)) and (flSl) . we find 

1 1 op 

544 = + + 2g) =[/ + / + 2 (-1 + /)] = -! + —, (50) 

2 2 r 
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where we have used Eqs. and (|^^ . In the weak- field approximation of GR, 
the component (744 is written in tlie perturbed form as 



544 



-1 - 20_ 



'N, 



(51) 



where 



GM 



m 
r 



(52) 



is the Newtonian potential, G, the gravitational constant, M, the mass of the 
source, and m = the geometrical mass. From Eqs. ([50)1 and ([5T|) . we 
conclude that 

F = m. (53) 

We conclude this discussion by stating that the metric of the vacuum static 
solution is completely specified by Eqs. gSl), (gH), (gS]) and This metric 

is given by 







V 











sin^ ( 









■m 

-1 + 2i 



\ 




4 + Hi 

m I 



(54) 



With this result, the Galilean trivial metric in spherical coordinates is obtained 
in the limit r 00 and it has a well defined relativistic limit. Eq. ([54]) will 
be called the Galilei-Schwarzschild solution. Notice that the event horizon, at 
r = 2m, is the same as in the relativistic context. However, this is of no concern 
for the applications discussed henceforth since the numerical values of the radii 
will be much larger than this critical value. 



3 Three classic tests of Galilean gravity 

In view of the Galilci-Birkhoff theorem, the Galilei-Schwarzschild metric, Eq. 
(|54|) is the only Schwarzschild-like metric consistent with Galilean covariance. 
The resulting line element, 

ds^ = -^ll-+r^d9^+r^smHdif^ + - {dx''f+- {dx'f+2 (-1 + -) dxUx^ 

^ r J- J- J- 

.... 

describes the invariant interval on the curved manifold in which the motion of 
particles and light rays will take place. In this section, we shall study these 
motions in order to test the Galilei-Schwarzschild solution. Along the lines of 
Ref. [8j, obtaining results from Galilean gravity that arc identical to results 
from GR means that these are not specifically relativistic in the Lorentz sense; 
in order the concept of curved manifold is sufficient to explain these effects 
without relying on the underlying Lorentz metric. 
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3.1 Planetary motion 

The motion of particles is described by the geodesic equations, 



ds ds 



For the Gahlei-Schwarzschild solution, these equations become 



1 



ds2 



1 - 



7.2 (jl _ 2m 

2m' 



r 



r sm 



^ - 1- 



+ -0, (56) 



2m 



m 
2^ 



2m 



dx^ 
ds 



dx^ 
ds 



ds"^ 



2drd0 
r ds ds 



dip 
ds 



3^ =0, 



d'^ip 2 dr dip 



dO dp 



-TT + — — r- + 2 cot — ^ = 0, 



ds 



2^4 



ds^ 



ds^ 



r ds ds 
1 



ds ds 



dr 



ds 
dr 



1 



ds 



dx^ 
ds 

dx^ 
ds 



dx^ 
ds 

dx^ 
ds 



0, 



= 0. 



(57) 
(58) 
(59) 

(60) 



We show that the motion lies in a plane as it happens in classical mechanics for 
a central force |12] . With an appropriate orientation of the axis, we can choose 
the initial conditions to be 



TT 

2' 



0, 



for some initial value of s. This choice implies that the motion of the test particle 
starts at the ecliptic plane with zero initial azimuthal velocity. It means also 



that Eq. l|57|) gives a zero initial elevation acceleration, 



= 0. Thus, at 







an infinitesimal proper instant later, 6'as = ^ and (^)^^ — 0, and similarly 
after another As. and so on. As a result, the motion is confined to the plane 



The geodesic equations, Eqs. (|56p to (pO)) . are then simplified to: 

2 / n \ / ; \ 2 



1 



ds2 



1 



dr\^ / 2m 
ds ) \ r 



dip 
ds 



1 m / 2m\ /dx^ da;5\ 



2r 



ds 



ds 



(61) 



d'^p 2 dr dip ^ 
ds^ r ds ds ' 



(62) 
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If we multiply Eq. ([62]) by , we find 

ds 



(65) 



where his a constant related to the conserved angular momentum of the particle. 
Similarly, by adding Eqs. ([55)1 and ([M)) . and by multiplying the result by 
fl - —V we are led to 



ds 



dx^ 
ds 



where k is another constant. By subtracting Eq. ([63]) from Eq. ([64]) . we find 

dx'^ dx^ 



(66) 



ds ds ' 
where b is constant. From Eqs. ([55)1 and (pT]) . we obtain 



ds 



and 



dx^ _ 1 
"ds" ~ 2 



fl _ 2m ^ 



(67) 



(68) 



(69) 



By substituting Eqs. 

d^r m 1 
d? 



and (inni) into Eq. ([HI), it follows that 



1 _ 2m^ \ 



dr 



1 



2m\ h"^ 1 



0. 



As in the classic Kepler problem, we can simplify the integration processes by 
considering r as a function of </? instead of s. If we change the variable r to 

1 

r 

it is possible to rewrite the last differential equation as 



d^u m / du 

d(p'^ {1 — 2mu) ydip 



+ (1 — 2mu) u 



2 /i2 (1 - 2mu) 



(70) 



The term proportional to (^^^ can be expressed in another form. For this, we 
use the constraint 

dx^ dx" 
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which leads to 



du , , 
— + (1 - 2mu) 
dip ' 



(72) 



By insertmg Eq. (j72p into Eq. jTOj), wc find the orbital equation, 



d'^u in 2 "^ 5^ 



0. 



(73) 



The first four terms are the usual ones obtained in the standard 4-dimensional 
GR. The extra term should provide non-relativistic corrections. 



Remark 1. Relativistic (Lorentzian) Limit. The projection of a 
Galilean theory to a Lorentzian theory can be performed by letting 

(x.«^.-H...,„.=).(x,-L._L). 

This is the case when, from Eq. ([3]), we have = x'^ and = 0. Therefore, 
in order to check consistency of our theory with the relativistic case, we should 
take = x^ , which implies = In this case, Eqs. ([68]) and ([69]) impose 
6 = and Eq. ((73)) reduces simply to: 

d^u m 2 

-r^ + u — —r — Amu = 0. 

This is the orbital equation obtained in standard GR for the Schwarzschild 
solution. So it is verified that the geodesic equation also has a good relativistic 
limit. 



Remark 2. Large distances. If the large distance case is considered 
then the terms proportional to v? in Eq. ([75]) should be neglected. Thus 

d^u m mb"^ ^ 

The first three terms are the usual terms obtained by Newtonian gravitation. 
The additional term, proportional to is open for interpretation. If agreement 
with measurements are to be obtained, then either h should be negligibly small 
(which would leave us with the usual GR result) or 6 should be included in a 
renormalized value of to. In Section [3?2l we shall use 5 << 1. 

3.1.1 The perihelion shift 

Let us rewrite Eq. ([73)l as follows: 

— J + u = + 3tou , (75) 
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with 

The orbital equation ([75|) is formally exactly as predicted by GR, the only 
difference being the redefinition h h. Then, we know before hand that our 
Galilean version of Schwarzschild solution will predict a perihelion shift for the 
orbit of the planets in the way of GR. 

The usual procedure is to obtain a solution of Eq. (j75p through an iterative 
procedure taking u ~ u^^^ +u^^^ [13| . The zero-order is the unperturbed solution 
of Eq. ((74|) . or Eq. ([75| with 3mu^ = 0. It is utilized as a source term for 

the differential equation of u^^\ i.e. we shall write 3mu^ = 3m (u'"') . The 
integration constants are the eccentricity of the orbit e and an arbitrary initial 
value ipo for the azimuthal angle. The constant e is related to the major axis 

a = rniax by 

where 

is the semi-latus rectum of the orbit. Considering orbits of small eccentricity 
(like the ones of Mercury), the terms of order and higher may be neglected 
and the solution of the orbital equation becomes 

TTl 

where 

Thus, while the planet rotates through an angle (p, the perihelion shifts by 
a fraction of the revolution angle given by 

AifQ ini? 3to 3m 



Lp L a (1 — e) 

where we have employed Eqs. ([78|) and ([77| . After a full revolution, ip ~ 2tt, 
we have 

„ m^ „ m „ GM 
Aipo = 67r^ = 67r- = Btt^— ^. 80 

3.1.2 Numerical analysis for the planet Mercury 

Here, the quantities of interest are expressed in terms of orbital parameters of 
the planet under consideration and the geometrical mass of the Sun. Let us 
consider the planet Mercury, with the following orbital data |15| : 

a = 57.91 X 10^ km, e = 0.2056. 
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The numerical value of the geometrical mass of the Sun m is calculated from 
the observational data: 



GMsun 132712 X 10^ W/s^ i .ysfi u ^o^^ 

m= 5 — = - -^j ^ = 1.4766 km, (81) 

c2 (299792.458) (km/s)2 

where GMsun = 132712 x 10^ kmVs^ [H]. 
Hence the value of L for Mercury is 

L = (1 - e^) a = (l - (0.2056)^) (57.91 x 10^) km = 55.462 x 10^ km 
and the perihelion shift, from Eq. (|80p . is 

A^o = 1-5974 TT x 10"^ 0.10351". 

This is an extremely small angle, but this is a secular effect which increases with 
the number of revolutions. The shift above is observed in a single Mercury- year; 
which corresponds to 0.241 Earth-years. So. the total shift per Earth-year is 

A^o 0.1035" „ 

AifiE = — ^ = = 0.42946". 

0.241 0.241 

If this effect is accumulated over 100 Earth-years, the total shift is 

A(^ = 100 Ai^B = 42.946" ~ 43". 

The conclusion is that our Galilean-Schwarzschild solution is as good as the 
GR Schwarzschild solution for a prediction of the perihelion shift of Mercury. 
The difference is that we can calculate the value for the constant h while in GR 
we obtain directly the value of h. Indeed, from the definition of the semi-latus 
rectum, given in Eq. (j78p . and the values of m and L, we obtain 

m 1.4766 o 

= — = = 2.6624 x 10"^. 

L 55.462 X 106 



Then, from (|76p . we find a constraint for the integration constant h and 6, 
11 



4 = 4(1 + ^) = 0.3756 xl0«, (82) 



so that 



1.6317 X lO"''. 



V0.3756 X IQS 

Note that our results arc derived by considering c as the speed of light (see 
calculation of geometrical mass of the Sun, Eq. (|5T|) ). 
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3.2 The deflection of light rays 



Light rays consist of massless test particles which travel at the speed of light. 
In special relativity, the photons move along the light-cone, following a null 
geodesic: ds^ = 0. They are essentially relativistic particles, and they will 
be treated accordingly here: we will keep ds^ = for photon traveling in our 
Galilean gravity background. So our framework is a curved manifold with local 
Galilean invariancc described by Galilci-Schwarzschild solution on which the 
ultra-relativistic particles will propagate. We feel justified in doing so because 
the photons are test particles and, by definition, test particles do not affect the 
geometry of the background spacetime (no matter their nature). 

By taking ds^ = in the Schwarzschild-like spacetime, Eq. ([55]) . and di- 
vinding the result by dcr^ (where a is an appropriate invariant length), the line 
element becomes 



1 f dr 
1- ^) 



-r^ sin^ ( 



dip 
da 



m 
r 



da 



dx^ \ 2 dx'^ dx^ ^ 
da I da da 



(83) 

This constraint replaces the one given by Eq. ((7T|) for massive particles. All the 
equations before Eq. (j7ip remain valid for the propagation of light, provided 
that we replace the invariant length s by a; that is. 



d r 
da^ 



1 



1- 



1 m 
2^ 



dr 
^da 

2m 



da 



2m 

r 

dx^ 
da 



dip 
da 



dip 
da 



dx^ 
da 

dx^ 
da 

where the initial conditions arc ( 



k 



(84) 
(85) 
(86) 

(87) 



9o = 7r/2 and (^)g = 0. These conditions imply 

and restrict our study to the plane 6 = 7r/2. Therefore, Eq. (|70p 

is still valid for the light rays. However, Eq. (jTS]) must be modified, because 
it was obtained using ds^ ^ 0, g^uU^u'^ = ^ 1, whereas we have here ds^ ~ 0, 
= 0. 

Let us rewrite the new constraint, Eq. (|83p . by substituting 6 — n/2, 
de/da = 0, together with Eqs. (gSl), dMl) and dST]): 



1 



2m 
r 



fc2 

2 



0. 
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Since 



we have 



drV 
dip] 



dr 
da 



dr dip 
dip da 



2m\ 



dr h 

dip r^ 



2 



7?2 



This equation can be written as a function of u = 1/r 

2 



du 
dip 



— + (1 - 2mu) 



lfc2 



(88) 



which differs only sHghtly from our previous constraint. Eq. ([71 
By substituting Eq. ([70)1 into Eq. one obtains 



dip 



iinu 



m 6^ 



(89) 



(90) 



If 6 ~ 0, we observe that Eq. (|89p reduces to the equation obtained in 
the standard 4-dimensional GR leading to the deflection of Ught rays. So it is 
verified that the equation for the path of light in the Galilean gravity has a good 
relativistic limit. 

In analogy with the previous subsection, we first set 3mu^ = 0, in order to 
get an approximate solution of Eq. (j89p by an iterative procedure, starting with 
a zero-order solution, 

(0) 1 / TO ^2 

-'^ = ^cos(^-^„)---, 

where ipo and R are integration constants. The interpretation of R becomes clear 
when we set ipo = 0, and introduce a Cartesian coordinate system x = r cos ip, 
y = r sin ip, with origin O at the center of the massive body, which is the source 
of the field. 

With & = 0, we find that Eq. (|90p reduces to i? = rcos(p = x. This is 
a straight line parallel to the y-axis, and R is the minimum distance between 
the light ray and the origin O. In this case, does not bend the straight 
trajectory of the photon and there is no defiection of light. However, this is 
not the best possible approximation, and there is also a 6 7^ to be taken into 
account. 



The first approximation to Eq. 



dip'^ 



,(1) 



TO 

7?Y 



3to 



becomes 



dip^ 



TO 6^ 
3to3 /b^ 



3to 



cos^ p 



— COSifi. 
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If we consider & << 1, then the factor of 5^ can be neglected, 

^ + ^ cos 6m--- cos 

A particular solution of this differential equation is 

fil ™ / 2 r, • 2 \ m b'^ m b^ 1 

" =;^(cos ^ + 2sin ^) - ;^Y-3TO-^--v?sm(p. 

Therefore, 



1 / 2 ^ • 2 \ 3to . \ 

- cos + (cos ^ + 2sin (^) - I 2 + — (^sin(^ 1 . (91) 



Notice that if 6 = this equation reduces to the expression derived in GR [13] ■ 
This also means that all possible modifications predicted by the Galilean version 
of the theory of gravitation for the deflection of light are encapsulated in the 
last term of Eq. jH]). 

If we multiply Eq. (|9ip by rR, we have 

R = r cos ifi + — (r cos^ (f + 2r sin^ <p) ~ tt i^ + (3mcp) r sin (p] . 
R ft,^ 2 



Then, in Cartesian coordinates, one has 



x^ R- 



2R\Jx'^ + + (^3m arctan — ^ ?/ 



The second term on the r.h.s gives the GR's deviation of the light ray from the 
straight line x = R. The last term is the contribution coming from the Galilean 
gravity. In the limit where y >> x, we obtain the asymptotic equation: 

in TfL bp" r / TT \ 

X'^R--{±2y) + —- 2i?(±y)+ 3m- y 



i? ' h'^ 2 



2 



because limQ,^+oo (arctan a) = 7r/2. Thus, the two possible values of x are 

2m m 6^ / Stt \ 

„ 2™ m b^ f „ 37r \ 
.-=i?-^y+^y(;-2i?+-mj,, 

and the deflection is described by the angle 

tano ~ . 

y 
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With tana = a + O {a^), the deflection angle is 



Am 

5=^ + {2mR)j^. (92) 

If wc take 6 = 0, wc consistently recover the GR's deflection angle, which 
acquires its maximum value when R is minimum, i.e. when the light ray coming 
from a distant source passes just outside the surface of the Sun. Thus, imposing 
that R is the radius of the Sun, R = 696,000 km and m = 1.4766 km [T3], we 
compute 

4?77. 

— = 8.4862 X 10"^ = 1.7504". (93) 
R 

The observations indicate that an average value for the deflection angle is about 
20% higher than this prediction: 

,5 = 2" = 9.6963 x 10"^. (94) 

Therefore, we observe that the extra degree of freedom arising from the Galilean 
gravity, 6, is such that it describes this difference between GR's prediction and 
observations. In other words, Eqs. ((92|) . ((93)) and ([94)1 can be used to constrain 
the value of the ratio b/h. Indeed, 

Now, in order to estimate the values of the parameters, consider Eq. ((82l) . 
which was deduced in the context of the shift of the Mercury's perihelion. By 
using the consistency of both constraints, Eqs. ([82]) and ([95]), we can determine 
b and h. This leads to 

1 / 5.8873 X 10-i3/i2\ ^^^^^ , o 
^ U + J = 0.3756 X 10^ 

from which we find ft., 

h = ^ = = 1.6317 X lO"'', 

^0.3756 X 10* - 

and 



b = V5.8873 X lO-13/i = 1.2520 x 10"^°. (96) 

This only gives a rough estimate of the parameters, because there is no formal 
connection with the parameters used in our previous treatment of Mercury's 
perihelion. 
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3.3 Gravitational redshift of spectral lines 



In this section, we examine the shift of the atomic spectral hnes in the presence 
of a gravitational field, also called the gravitational redshift. 

From the Galilei-Schwarzschild line element, Eq. (|55p . we define 



dr' 



--ds^ 
2 



as the time-interval between two events with vanishing spatial separation, dr = 
dO = dip = 0. The minus sign comes from our choice of the signature of the met- 
ric, and the factor 1/2 is chosen in order to allow agreement with the relativistic 
analog when b = 0. 

The time- interval is related to the "coordinate time" differential, dx'^, and 
the additional Galilean differential coordinate, dx^, by 



_ 2dT^ = ^ {dx^) 



{dj 



2-1 



dx'^dx^ 



(97) 



where the differentials are constrained by Eqs. ([55]) and 



dx^ 
d7 



(1 



dx^ 
77 



so that 



dx^ = dx^ - bdr. 



(98) 

Note that these last three equations are very general; they are valid for massive 
particles, but they have exactly the same form as for massless particles, such as 
photons. 

With Eq. ([55]) . the time-interval in Eq. ([W]) can be expressed as 



dr^={l- ^) (dx^) 



Notice that if 6 = 0, we have 



dT 



— dr" 

r 2 



1 - 

r 



1 



2 TO 

r 



bdx dr. 



(99) 



(100) 



which is the expected Schwarzschild solution of GR (see Ref. [TS], Eq. (4.92)). 
The positive sign follows from the natural assumption that the proper time r 
should increase with the time coordinate x^. With 6^0, Eq. ([99|) leads to 



%x^Y = 0. 



This is a second-order equation for dr. By solving for dr = dr (rfa;^), we find 



dT 



{l-^)b dx^ ±^{l-2j^fb2 (^^4)2 ^ 4 (1 ^ (1 _ 2j^^ 

2(1 + Tt) 
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In order that this expression agrees with Eq. (jlOOp for 6 = 0, we must choose 
the plus sign in the r.h.s., which leads to 



dr 



1 + ^- 



I _ 2™. 



2m 



(101) 



This distinction between proper time and the time coordinate gives rise to a 
difference between the proper frequency and the coordinate frequency of a peri- 
odic phenomenon in the curved spacetime, like the emission of electromagnetic 
radiation by an atom. Consider the propagation of electromagnetic waves in the 
limit of geometrical optics. Then, the electromagnetic field can be written as 
f = a exp {iip) , with a = a (x, t) the wave amplitude, and the phase ip = {x, t) 
an eikonal function. Then the frequency of the wave can be expressed as the 
derivative of ip with respect to the time, and one has a coordinate frequency, 
ojo = and a proper frequency, oj = 

We call LUi the proper frequency of the wave emitted by an atom at a point 
Pi. At another point, P2, the observed proper frequency will be different, say 
L02, once the gravitational field is not the same. They are related so that 



UJ2 



LUl I dip 
dr 



dt \d 



■7)1 ^o(arji 



with ijjQ constant and where is the time coordinate t. From Eq. (jlOip . we 
find 



dx^ 



1 -X- 2m 62 
+ n 4 



1 



1 I & 

^ 4 2 



1 _ 2m /x _ 2m 
ri y ri 



and similarly for {dx'^ / dr) ^. Then, the ratio of frequencies is 



UJ2 



1 -L 2m 
^ r2 4 



1 -I- 2m 6^ 
^ ri 4 



1 + ^ _ & 
^4 2 



Y 2m 

ri 



4 



2 _ 2m 
ri 



2 _ 2m 
r2 



(102) 



This exact result can be approximated by taking into account the fact that 
regular estimates assume ri << r2, and that b is small (see estimate in the 
previous section, Eq. ([96]) ). For instance, let us assume that the radiation 
observed on Earth, at r = r2 = 149.6 x 10^ km, was emitted at the surface of 
the Sun, at r = ri = i? = 696,000 km. Then the approximation ri << r2 is 
certainly valid. Moreover, the quantity m/ri is very small, since m = 1.4766 
km, and the functions in Eq. (|102p can be approximated accordingly. Therefore, 
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for m/ri « 1 and m/r2 << 1, the ratio of frequencies is 



1 + ^ - ^ 



^ - ^ ( 1 - 

4 2 I 7-2 



_^ TO 
J"! 



TO 

1 + — 

?'2 



If we assume that m/r2 — 0, then this ratio simphfics to 



UJ2 



1 - 



2to fe2 

ri 4 



1 + £1 _ b 

^ 4 2 



rn 



4 



1 - 



Furthermore, if h is also assumed to be very smaU, 6 << 1, this equation becomes 



UJ2 



TO 

ri 



(103) 



so that 



LU2 ~ ^1 



m 



1 - 



When 5 = 0, this result agrees with GR. 

From Eq. (|103p . the frequency iU2 measured on Earth is smaller than the 
proper frequency, uji , and we obtain that the spectral line produced on the Sun 
are displaced towards the red compared with the same line produced on Earth. 
With the values of to, R and b utilized in the previous section. 



CO 



1.4766 
696000 



1 



1.2520 X 10 



-10 



-2.1216 X 10" 



This is too small to be measured and the correction due to b is completely 
negligible. One should keep in mind that this is simply an estimate based 
on the parameters obtained in Eqs. ([5^ and Accordingly, everything 

seems consistent with GR, which shows that the gravitational redshift is not 
specifically relativistic in the sense of local Lorentz invariance. 

Let us emphasize that this effect can be verified by observing massive white 
dwarfs, or with a Mossbauer experiment designed to detect the deviation in the 
wavelength of the radiation emitted in the presence of the gravitational field of 
the Earth. 



4 Concluding remarks 

In this paper, we have constructed an exact Schwarzschild-like external solution 
of the gravity field equations based on local Galilean invariance. In a flat space- 
time, Galilean covariance consists in writing dynamical equations in a (4 + 1) 
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Minkowski manifold with light-cone coordinates, thereby providing a unifying 
treatment of relativistic and non-relativistic field theories. This is extended 
here to curved spacetimes. We note that our exact solution comprises an event 
horizon which turns out to be much smaller that the distances involved in our 
numerical examples. Despite being exact, the Galilei-Schwarzschild solution is 
said to be effective in the sense that it would not be suitable to describe near- 
blackhole situations, in the very strong gravity regime. We examined three GR 
classic tests: the perihelion shift of Mercury, the deflection of light by the Sun, 
and the gravitational redshift of atomic spectral lines. The Galilean theory gives 
results similar to the GR predictions. This means that the results remain the 
same, whether covariant theory of gravity is locally Galilean invariant or locally 
Lorentz invariant. 

Possible avenues of research include: (i) the search for non-static solution 
(i.e. a;^-dependent) like Kerr-Galilei metric or Friedmann-Galilei metric. An 
issue with this last proposal is to determine which coordinate should be singled 
out as the time coordinate in order to define a cosmic time, for instance; (ii) the 
analysis of the Killing vectors for Galilei-Schwarzschild metric and the study 
of its invariants; (iii) the investigation of MOND models [16] in the context of 
Galilean gravity, and applications to the description of the cold dark matter and 
the rotation curve of spiral galaxies; (iv) the relation between our parameter 
b and the parameters appearing in the different post-Newtonian expansions; 
(v) the construction of the gauge theory of gravity based on local Galilean 
invariance. 
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